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Consequences of the constraint of SU(3) colour singletness 
of quark-gluon plasma are studied. This restriction increases 
the free energy barrier for the formation of hadronic bubble in 
supercooled phase and influences significantly the dynamics 
of the initial stage of quark-hadron phase transition. It also 
introduces terms dependent on the volume occupied by the 
plasma in the energy density and the pressure. These mod- 
ifactions vanish in the limit of an infinite volume. The last 
stage of the hadronization of the QGP likely to be formed 
in relativistic heavy ion collisions is necessarily characterized 
by a decreasing volume containing the quark matter, and 
thus these corrections become important. The nucleation of 
plasma droplets at AGS energies is also seen to be strongly af- 
fected by the requirement of colour singletness, and the choice 
of prefactor. 

PACS numbers: 24.85. +p, 12.38.Mh, 25.75.+r, 64.60. Qb 



I. INTRODUCTION 

The success of the quark model, the quantum chromo- 
dynamics (QCD), and the non-observabiUty of the free 
partons (q, q, g) has entailed the concept of confinement. 
QCD, the theory of strong interactions, is not perturba- 
tive at large distances. Thus, the confinement itself can 
not be treated perturbatively. There are reasons to be- 
lieve that the confinement of partons inside hadrons may 
not survive collisions between heavy nuclei at relativistic 
energies. In such collisions, the two nuclei masquerading 
as clouds of space and time like partons pass through each 
other, leaving behind a high density plasma of quarks, 
antiquarks, and gluons (QGP) [Q in their wake, in the 
region between the two receding fronts of the leading par- 
ticles. This plasma expands and cools, the energy density 
becomes low enough and a phase transition to a hadron 
gas takes place around the critical temperature, Tq. A 
dynamic treatment of this phase transition is a problem 
of considerable interest. 

An understanding of the QCD phase transition re- 
quires a knowledge of the equation of state as well as 
the kinetics of phase transition. If QCD has a first order 
phase transition, it may proceed with a supercooling of 
the QGP followed by a nucleation and growth of hadronic 



bubbles |g,^ , releasing the latent heat as the phase tran- 
sition progresses. In a superheated hadronic matter, on 
the other hand, nucleation of a QGP droplet may also 
proceed similarily. 

For a first order phase transition the rate for hadronic 
bubble/plasma droplet nucleation can be estimated in 
the frame work of homogeneous nucleation theory M 



/ = /ocxp(-AFVT) 



(1) 



where Jq, which has the dimensions of 1/fm^-is called the 
prefactor, T is the temperature, and AF* is the change 
in free energy accompanying the formation of a critical 
size hadronic bubble/plasma droplet. These dimensional 
arguments were used in a large number of studies in 
the past to replace /q with T* or Tq. This unsatisfac- 
tory state of affairs was corrected recently by Kapusta 
and Csernai |g]. They computed the dynamical prefac- 
tor in a course-grained effective field theory approxima- 
tion to QCD. This dynamical factor influences the growth 
rate and statistical fluctuations and also accounts for the 
available phase space. 

Csernai et al. [pP] have also used a nucleation rate 
equation with this realistic dynamical prefactor to study 
the time evolution of expanding QGP as it converts to 
hadronic matter. They noted a substantial deviation 
from an idealized Maxwell construction that has often 
been employed as a model of hadronization l6l. Obvi- 
ously, such an idealized phase transition assumes a QCD 
nucleation rate which is much larger than the rate of ex- 
pansion. This is not necessarily true. 

In all such studies, QGP is generally described as an 
ideal gas of quarks, antiquarks and gluons, essentially 
described by the Stefan-Boltzmann law. Lattice calcula- 
tions have provided ample evidence that even at fairly 
high temperatures, colour singlet objects like multi-quark 
cluster (qq, qqq, qqq, ■ ■ •) propagate in the plasma. One 
may account for this 'interaction' by requiring that all 
physical states be colour singlet with respect to the SU(3) 
colour gauge group ^-Q. 

It has recently been shown |jl3] that restricting the 
quark partition function to be colour singlet of SU(3) 
colour gauge group amounts to reordering the thermo- 
dynamic potential in terms of the colourless multi-quark 
modes {qq, qqq, qqq, ■ ■ ■) at any given temperature. Un- 
der a suitable confining mechanism, these could evolve 



into colour singlet hadrons/baryons at low temperatures. 
This is also in accord with the "preconfinement" prop- 
erty of QCD noted by Amati and Veneziano |1^ quite 
sometime ago where the cascading and fragmenting par- 
tons produced in hadronic collisions rearrange themselves 
into colour singlet clusters which ultimately evolve into 
hadrons ||lj,|lj| ■ These considerations convince us that it 
is important to incorporate the dynamic requirement of 
colour singletness of the quark-matter which "tunnels" 
into hadronic matter phase space |15|| . 

In the present work we study the consequences of 
the incorporation of colour singlet equation of state for 
plasma on the dynamics of quark-hadron phase transi- 
tion. 



II. EQUATION OF STATE 

A. Colour Singlet Equation of State for a 
Quark-Gluon Plasma 



We note that the colour singletness introduces correc- 
tions to the normally assumed expresssions for the energy 
density and the pressure which vanish for an infinite vol- 
ume. Now consider the hadronic phase arising as a result 
of first order quark-gluon/hadron transition through the 
nucleation of hadronic bubble in the bulk QGP. With 
the passage of time, more and more of the quark matter 
will get converted to hadronic matter. The later stage of 
the process of hadronization will be characterized by a 
decreasing volume occupied by the quark matter. If we 
believe the above equations of state, the volume occupied 
by the plasma can not be vanishingly small. The colour 
singletness will again have important consequences for 
dynamics of the phase transition during the later stage 
of hadronization. 

Here, we would also like to point out that the colour 
singletness has important bearing on the nucleation of a 
hadronic bubble in the plasma which we will sec in sec. 
III. 



Consider a quark-giuon plasma consisting of 'u' and 'd' 
quarks, and gluons. The grand canonical partition func- 
tion p-lllll subject to colour singletness can be written 
as 



Z{/3,Vg)^TY (Ve-^'" 



(2) 



where f3 — 1/T is the inverse temperature, Vq is the 
volume, H is the Hamiltonian of the physical system, 
and V is the colour projection operator. For a baryon 
free plasma, this can be simplified after a considerable 
amount of group theoretic algebraic manipulations ppT| 
to give, 
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where a^ = 37tt^/90. Now the free energy of the baryon- 
free colourless quark-gluon gas is obtained as. 



F, = 



T liiZ{T,Vq) +BV, 



1 ■ 



(4) 



One may now write for the energy density 
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The pressure of the above quark-gluon system is given 



as 
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dF 
dV 
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B. Equation of State for Hadron Gas 

We model the hadronic phase as a gas of massless pi- 
ous. The energy density and the pressure of such a sys- 
tem can be written as 



and 



where ah — 7r^/30. 



e/i = iahT 



Ph - ahT^ , 



(7) 
(8) 



III. SUPERCOOLING AND NUCLEATION 

Nucleation in a pure phase like QGP proceeds via cre- 
ation of a hadronic bubble due to statistical fluctuations 
in a supercooled plasma. The bubble is made up of hot 
pion gas and is surrounded by colour singlet plasma of 
volume, Vq — {V — Vb)^ where V is initial volume of 
plasma. Vf, — AttRi, /3 represents an excluded volume 
corresponding to hadronic bubble. One can think of curv- 
ing out a colour singlet piece of palsma and replacing it 
with a hadronic bubble. The fields in plasma obey the 
bag boundary conditions, staying outside the hadronic 
bubble |16| , p7t . If the radius of the bubble is Rb, the 
change in free energy can be written pT| , p^ within the 
bag model as. 



AF = Tin (ttVs) -f 4Tln ( -VbT^ ] + aqVbT 



-{B + Ph)Vb + AirRb^a 



(9) 



where P^ is the pressure of the hadron gas given in 
Eq.(@) and a is the surface free energy of the quark- 
gluon/hadron interface. The first two terms in Eq.(|^) 



are due to SU(3) colour singlet restriction. They increase 
the barrier for AF required to form hadronic bubble in 
plasma (Fig.l). We shall see later that it has an impor- 
tant effect during the initial stage of QCD phase transi- 
tion. 

Recall that one can derive the critical radius {Rb*) of 
hadronic bubble by minimizing the change in free energy, 
A_F, with respect to i?(,. If they are too small {R < Rb*), 
they will shrink and vanish. If they are large {R > Rb*), 
they will grow. Now, minimizing AF with respect to Rb, 
one gets 



12T 

—^ - A7r{B - aqhT^)Rb 

Kb 
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STTcrRb* = , 



(10) 



which will yield the critical radius of hadronic bubble 
Rb*. We have further defined, a^/j = aq~ a^. If the first 
term in Eq.(nQ), which has its origin in the requirement 
of colour singletness is neglected, one obtains the criti- 
cal radius, Rb* — 2a /{B — UghT*), for a nonsingiet case. 
Making a substitution a; = l/Rb*, Eq.(p^ can be written 
as 







(11) 



where a = 2Tra/3T and b = tt{B - aqhT*)/3T. The 
physical solution of Eq.(ll) gives the critical radius of 
hadronic bubble for the colour singlet case, as. 



Rb* = 'Sz/ (Sz^ - 



6/2+ VaV27 + P74 



(12) 



Now, the change in free energy for the creation of a 
hadronic bubble having the critical radius Rb* is 



AF, = AF\r=r, 



(13) 



The rate for the nuclcation of hadronic phase out of 
the plasma phase is usually estimated by Eq.(|i|) with 
prefactor, Iq. As remarked earlier the prefactor, /q, has 
been calculated by Csernai and Kapusta |^] in an effective 
field theory approximation to QCD as 



/ = 16 f a y/^ TVgRb* 



Stt \3TJ ^^{Aw) 



(14) 



Here, rjq and £,q are, respectively, the shear viscosity and 
correlation length in the plasma phase, and Aw is the 
difference in the enthalpy densities of the two phases. 
We use the same parameter set as used in Refs. [^||, 
e.g., 51/4 = 235 MeV, ^g = 0.7 fm, 77, = UAT^ and 
Tc — 169 MeV. Next we closely follow the arguments of 
Rcfs. Ig,^ to obtain the dynamics of the phase transition. 
Once the nucleation rate is known, one can calculate 
the (volume) fraction of space h{T) converted to hadronic 
gas at a proper time r. This proper time is measured in 
the local comoving frame of an expanding system. For 
this purpose one needs a kinetic equation which involves 



/ as the source. If the system cools to Tq at a time tc, 
then at some later time r the fraction of the space which 
has converted to hadronic gas is |g,^ 



h{T) = / dr7(r(T'))[l-/i(r')]V(T',T) 



(15) 



where, V{t',t) is the volume of a bubble at a time t 
which has nucleated at an earlier time r'. This also takes 
into account bubble growth. How rapidly does the bub- 
ble grow after nucleation? Usually a critical size bub- 
ble is metastable and will not grow without a perturba- 
tion. Pantano and Miller |1§[| have numerically computed 
the growth of bubbles using a relativistic hydrodynamics. 
The asymptotic radial growth velocity was found to be 
consistent with the growth law 



viT) =. t;o (1 - T/Tc)'/' 



(16) 



where vq is a model-dependent constant. We shall use 
vo — 3c as it has been argued |g,|| that the Eq.(|l6|) is 
intended to be applied only as long as T > 2Tc/3 so that 
the growth velocity stays below the speed of sound of a 
massless gas, c/\/3- Thus the growth of the bubble can 
be approximated as, 



V{T',r) 



Rb*{T{T')) 



At"v{T{t")) 



(17) 



Now one needs a dynamical equation which couples 
the time evolution of the temperature to the fraction of 
space converted to hadronic gas. For this purpose we use 
both Bjorken longitudinal hydrodynamics and Cooper- 
Frye-Schonberg spherical hydrodynamics as Refs. [^,||. 
In Bjorken model [M the time evolution of energy den- 
sity e is given as 



de 
d7 



w 

T 



whereas in Cooper-Frye-Schonberg model |E0[ it is 



de 
d7 



Zw 

T 



(18) 



(19) 



In Eq.(pj|) the factor of 3 appears because of the spa- 
tial expansion along three dimensions rather than one 
dimension. In order for an easy comparison we con- 
sider initial conditions as in Refs. |0,p[ which are likely 
to be achieved in collisions involving two gold nuclei at 
RHIC energies. Thus for longitudinal expansion we take 
Ti = 3/8 fm/c and Ti = 2Tc- The temperature decreases 
as T{t) oc r"^/-^ until the time tc = 3 fm/c when the 
temperature becomes Tc- In the case of spherical ex- 
pansion we consider Ti = 2Tc and r^ = v2^nuci and 
the temperature decreases like T(t) oc r"'^ until the time 
Tc — 18 fm/c corresponding to half density radius of a 
gold nucleus (For details see Refs. ||]). The Eqs.(|l|) 



and (19) are essentially the statement of energy conserva- 
tion in the respective picture which assume kinetic equi- 
librium but not phase equilibrium. The energy density 
can be written I2IJ3I as 



3(T) = h{T)e^,{T) + [l-h{T)]e,(T) 



(20) 



Here, eh and Cg are the the energy densities in hadronic 
and plasma phases at temperature T, and similarly for 



IV. DYNAMICS OF PHASE TRANSITION: 

LONGITUDINAL AND SPHERICAL 

EXPANSION 
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In Fig.l the free energy difference, AF in Eq.(g) as 
a function of hadronic bubble radius for a fixed bag 
constant {B^^'^ = 235 MeV), surface tension {a = 50 
MeV/fm^) and temperature (T = 160 MeV) is displayed. 
Results are given with and without colour singlet require- 
ment of the system. One finds a significant increase in 
AF if the constraint of colour singletness is imposed. 
We will see below that it has important consequences on 
the time evolution of expanding QGP as it converts to 
hadronic matter through nucleation of hadronic bubble. 
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FIG. 1. The free energy difference AF(R) for creation of 
a hadronic bubble in quark-gluon plasma. 



In Fig. 2 we show the variation of the temperature with 
proper time as the matter undergoes longitudinal expan- 
sion with the initial conditions given in sec. 3. We have 
also given the results for an adiabatic phase transition for 
a comparision. The value of a used here is 50 MeV/irn^. 
The general features are similar to those of Refs. P,p|. 



FIG. 2. The temperature as function of proper time for a 
hadronizing quark-gluon plasma in a central high-energy nu- 
cleus-nucleus collision for matter undergoing longitudinal ex- 
pansion. The initial conditions correspond to QGP formation 
at BNL - RHIC energies. 



The matter continues to cool below Tq until nucleation 
of hadronic bubble sets in. For colour singlet case the de- 
gree of supercooling is about 30%, i.e., 12.5% more than 
that for the nonsinglet case. Once nucleation and growth 
of bubble start, the system reheats near Tc due to re- 
lease of latent as the phase transition progresses. When 
temperature approaches ~0.95rc nucleation of further 
bubble formation ceases off and the transition proceeds 
because of growth of previously nucleated bubbles. How- 
ever, the system can not reheat upto Tc because bubble 
growth (Eq.(|l7|)) reduces to zero as Tc is approached. 
Recall that nucleation in colour singlet case will be de- 
layed due to increase in height of AF which in turn also 
slows down the phase transition. This results in 10% ex- 
tra entropy generation in the processs as compared to 
colour nonsinglet case. 

Fig. 3 shows the variation of critical radius as a func- 
tion of proper time for hadronic bubble undergoing lon- 
gitudinal expansion with and without colour singlet re- 
quirement. It clearly indicates that the early stage of 
nucleation is characterized by a much larger critical size 
of nucleated hadronic bubble (~ 1 fm) for colour sin- 
glet case as compared to ~ 0.5 fm for colour nonsinglet 
case. At later times (> 12 fm/c) the critical size of nu- 
cleated hadronic bubble remains lower until the phase 
transition is completed if the colour singlet requirement 
is implemented in the system. In these initial studies we 
have not included the modification due to bubble fusion 
which should be of interest. 
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FIG. 3. Radius of critical hadronic bubbles as a function 
of time, in a hadronizing quark-gluon plasma. 
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FIG. 4. The volume-fraction of space occupied by the 
hadronic matter as a function of time, in a hadronizing 
quark-gluon plasma. 



In Fig. 4 we show the volume fraction converted to 
hadronic matter undergoing longitudinal expansion as 
a function of proper time. The delay in completion of 
QCD phase transition due to colour singletness is seen 
clearly. As the phase transition progresses via nucle- 
ation of hadronic bubble, the available space will be 
progressively occupied by hadronic matter. We would 
like to make an amusing observation here. If we believe 
the equations of state for the QGP as discussed earlier 
(Eqs.(||) and (^)), the volume occupied by plasma can 
not be vanishingly small. This necessarily implies a rem- 
nant of quark matter when the process of hadronization 
is over. We find that a fraction of quark matter of mass 



~ (5-10) GeV having volume ^ 10-15 fm^ remains un- 
converted at the end of the hadronization in the central 
region, per unit rapidity. One can have interesting spec- 
ulations about such a remnant. 
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FIG. 5. Same as Fig. 2 with various a values for colour 
singlet case. 



In Fig. 5 we have attempted to study the dynamics 
of phase transititon as a function of quark/hadron in- 
terface tension a. Recall [pi that the dynamics depends 
sensitively on a when we do not impose the restriction of 
colour singletness. With imposition of colour singletness 
large variation in a leaves the dynamics fairly unchanged. 
The reason for this is not too difficult to identify. Note 
that Eq.(|o|) tells us that even if a = 0, there is pos- 
sibility of nucleation of hadronic bubble with a critical 

1 /^ 
radius, Rb* = [3T/7r{B — aqh)] ■ However as the nu- 
cleation rate, /, itself depends on a via its prefactor, Iq 
and / == if cr = 0. Thus we have kept ct = 50 MeV/fm^ 
as / involves a more realistic prefactor characterized by 
a. 

How will these findings differ for a more realistic (3+1) 
dimensional expansion of plasma? Normally, the plasma 
is expected to expand mostly in longitudinal direction 
initially. After a time r ~ R/c^ where R is transverse 
radius and c^ is the speed of sound, the system is likely 
to expand in transverse direction as well. We, like the 
authors of Refs. M, take the other extreme and look on 
spherical expansion with the same initial conditions as in 
Ref. @. 
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FIG. 6. Same as Fig. 2 for matter undergoing spherical 
expansion. 
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FIG. 7. Same as Fig. 3 for matter undergoing spherical 
expansion. 



Fig. 6 is similar to Fig. 2 with the exception that it 
considers spherical (3-dimensional) expansion instead of 
longitudinal (l-dimcnsional) expansion. The degree of 
supercooling is almost of the order of longitudinal one 
(30%) with the prime difference that hadronization is 
faster in 3-dimensional case. This implies that the sys- 
tem spends less time in the neighbourhood of Tc than 
otherwise. Fig. 7 shows the variation of critical radius 
of hadronic bubbles undergoing spherical expansion as 
a function of proper time with and without colour sin- 
gletness. It is clear that the critical radius of nucleated 
hadronic bubbles for spherical expansion is smaller com- 
pared to longitidunal expansion. In Fig. 8 we show the 
variation of volume fraction converted to hadronic mat- 



ter for spherical expansion as a function of proper time. 
This is also clear from this figure that the phase transi- 
tion is faster than the longitudinal one. For the shake 
of completeness we also present Fig. 9 for spherical ex- 
pansion which shows the variation of temperature as a 
function of proper time for different quark/hadron inter- 
face values. 
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FIG. 8. Same as Fig. 4 for matter undergoing spherical 
expansion. 
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FIG. 9. Same as Fig. 5 for matter undergoing spherical 
expansion. 



V. NUCLEATION RATE OF DROPLETS OF 

QUARK-GLUON PLASMA IN HOT HADRON 

GAS 

The creation of QGP at AGS energies may pro- 
ceed through the nucleation of a plasma droplet in 
hot hadronic gas. Here the fields in plasma obey the 
bag boundary conditions, staying inside the plasma 
droplet ||l6|. If the radius of the droplet is Rq, the nucle- 
ation process is, generally, activated by the change in free 
energy which can be written within the bag model [|liyi6| 



AF = Tin [nVs] + 4Tln I -VqT^ \ - aqVqT* 

+ {B + Ph) Vq + AuRq^a . (21) 

Here, Vq is volume of the plasma droplet formed in a su- 
perheated hadron gas. Usually, the nucleation rate of a 
plasma droplet from superheated hadronic phase is esti- 
mated from Eq. (|l|) . Now, the prefactor [^ is given by 



/o - ^^0 



(22) 



The dynamical prefactor k, determines the exponential 
growth rate of critical size droplets, and is given by pi] 



2cr 



{tlwYRq*' 



AT + 2 ( -r; + C 



(23) 



where A is the thermal conductivity and ry and C are the 
viscosities of the hadronic phase. The bulk viscosity C, 
is very small compared to shear viscosity rj and can be 
neglected. For these dissipative coefBcients we use the 
parameterization of Danielewicz [E2|. Rq* is the critical 
radius of a nucleated plasma droplet. This can be ob- 



tained by minimizing AT, given in Eq.(21), with respect 
io Rq. To a first approximation the statistical prefac- 
tor RpTI is given by 



n. 



= J_(Z.\"^ (^ 



3V3VtJ \0i 



(24) 



For our purpose we use following set of parameter a — 
50 MeV/fm2, ^^ = 0.7 fm and B^l^ = 200 MeV which 
gives Tc '--^ 170 MeV. We shall also give results for the 
prefactor ~ T^ which is used often in such studies. 

As a first step, in Fig. 10 we plot the variation of AT as 
a function of droplet radius with and without colour sin- 
glet restriction for three different temperatures (160, 180, 
240 MeV). Due to colour singlet restriction, the height 
of AT is enhanced significantly in each case. Fig. 11 
shows the variation of critical radius of nucleated plasma 
droplet as a function of temperature. We see that the 
imposition of the colour singletness increases the critical 
radius by factor of ~ 3-5 over the range of temperatures 
that we have considered here. Knowing that creating a 



bigger size bubble by statistical fluctuation is consider- 
ably less probable, it is not surprising that this should 
lead to considerable suppression of nucleation rate when 
the degree of superheating is small. 



2500 



2000 



500 



-500 



- — No Colour Singlet 


' 


Colour Singlet 

B''''' = aOO MeV ^''~ 


_T=160 MeV 


^, 


17 = 50 MeV/fm^ / 


\ 


/ 180 


\ ■ 


240^ \ 


Y 


., - - xpa \ ~~ -^ \ 

\ 240 s 1 ■- 

\ K 


\^, , 



0.0 0.5 1.0 1.5 2.0 2.5 

R (fm) 



FIG. 10. The free energy diflerence AF{R) for creation of 
a plasma droplet in hot hadronic matter. 
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FIG. 11. Critical radius of nucleated plasma droplets as a 
function of temperature in a superheated hadronic matter. 



Fig. 12 shows the variation of nucleation rate of plasma 
droplets in superheated hadron gas as a function of tem- 
perature. It is seen that the nucleation rate with the pref- 
actor T^ (dashed lines) is suppressed considerably when 
the restriction of colour singletness is imposed. Similar 
finding has been reported recently in the literature p3] 



with a somewhat different expression for AF. The results 
(sohd Unes) with the more reaUstic prefactor (Eq.(|2^)) 
are richer in detail. Firstly we notice that without the 
restriction of colour singletness the traditional prefactor 
considerably overestimates the rate of nucleation. With 
prefactor given in Eq.(22), we see that switching on the 
requirement of colour singletness lowers the rate of nu- 
cleation at smaller T . However, in a surprising finding, 
we note that while the no-colour-singlet rate decreases 
with increase of T, that with colour singlet restriction 
increases and ultimately becomes larger at really high 
T. This interesting behaviour is seen to emerge from 
the structure of the prefactor (Eq.(p2[)) which in fact is 
proportional to Rq* , and Rq* is seen to decrease with 
increase in T (Fig. 11). We add here that Kapusta and 
Vischer |24] have recently proposed a more efhcient mech- 
anism for nucleation of plasma droplet at AGS energies 
which envisages seeding of the plasma via collision of two 
very energetic nucleons in the hot and dense hadronic 
matter. 
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FIG. 12. The nucleation rate of plasma droplets as a func- 
tion of temperature, in a superheated hadronic matter. The 
solid lines are with dynamical prefactor of Ref. fej whereas 
those with dashed lines are with T^ as prefactor. 



VI. SUMMARY 

Before summarizing let us briefly examine some of the 
inputs in the present work. While writing the expressions 
for the gain/loss in free energy (AF) we have included 
only the volume and the surface terms. This should be 
valid when the size of the plasma droplets or hadronic 
bubbles is large. The prefacor we have used is valid for 
the situation when the size of the bubbles/droplets is 
larger than the correlation length of the system (^ = 0.7 
fm) . We have seen that these bubbles/droplets have radii 



larger than a fermi, and thus both these conditions are 
reasonably satisfies. 

Recall that the surface energy coefficient a is zero in 
the MIT bag model for massless quarks |2^. Quarks 
traversing a hot medium, do acquire a thermal mass. 
Finite temperature lattice QCD |g6| and a pure SU(N) 
gauge theory [g^ yield a value in the range a « 20-70 
MeV/fm^. We have thus mostly used cr w 50 MeV/fm^ 
and the effect of varying a is also studied. We have not 
come across any estimate of the so-called curvature term 
in this case. We have verified, however, if we add a term 
for curvature energy in the expression of AF for the cre- 
ation of hadronic bubble in the plasma as in Ref. p6J, 
then the supercooling of the plasma is completely el- 
eminated, if we ignore the colour singletness. We have 
also checked that supercooling of plasma is possible if 
the colour singlet restriction is imposed along with cur- 
vature contribution. However, for plasma droplet in hot 
hadronic gas the degree of superheating will be enhanced 
than otherwise. 

In brief, the theory, proposed recently ||] , to describe 
the dynamics of hadronization has been generalized to 
study the consequences of the requirement of colour sin- 
gletness of QGP which may be produced in relativistic 
heavy ion collisions. It is shown that hadronization of 
longitudinally and spherically expanding plasma may be 
slowed down due to this requirement. While there is 
no production of entropy for an adiabatic phase transi- 
tion, the entropy increases by about 30% in the treat- 
ment of Ref. yj. The requirement of colour-singletness 
introduced in the present work enhances this entropy by 
an additional 10% and also increases the degree of su- 
percooling. We note an interesting possibility; a small 
fraction of QGP may not hadronize at all. 

We also find that the nucleation of a plasma droplet in 
a superheated hadronic matter is suppressed at low tem- 
perature (T < 220 MeV), but it is enhanced at higher 
temperature (T > 220 MeV) when the colour singlet re- 
striction of QGP is accounted for. The decrease at T '--^ 
170 MeV is found to be by four orders of magnutude 
while the enhancement at 250 MeV is by one order of 
magnitude. 

It would be of interest to consider a generalization of 
this study to the case of non-equilibrated plasma which 
hadronizes into a hadronic matter having a richer equa- 
tion of state. A (rather) crude result may be obtained 
by keeping all other parameters like Tc-, cr, r/, etc. fixed 
to their present values. This we feel, may not be quite 
justified. 



We gratefully acknowledge useful discussions with Joe 
Kapusta during this work. 
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